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ABSTRACT 



A test facility to evaluate the effect of condensate 
inundation on heat transfer within a horizontal tube bundle 
was designed, constructed and validated. Five 15.9 mm (5/8 
in.) nominal outside diameter, smooth stainless steel tubes 
were utilized in a vertical row. They were located in an 
equilateral triangular array with a spacing to diameter ratio 
of 1.5. 

Heat transfer performance was determined for each tube 
in the bundle. Data was taken by condensing steam at about 
21 kPa (3 psia) on the outside of each tube. Each tube was 
cooled by water on the inside at velocities of 0.78 to 7.0 
m/sec (2.56 to 23 ft/sec). The overall heat transfer coeffi- 
cient was determined directly from experimental data. The 
inside and outside heat transfer coefficients were determined 
using the Wilson plot technique. 

Observation of condensate flow showed lateral droplet motion 
along the tube in portions of the condenser as well as side 
drainage, particularly over the first three tubes. Outside 
heat transfer coefficients were lower than expected when com- 
pared to Nusselt theory, possibly due to the effects of 
secondary vapor flow and/or non-^-condensable gases. Recommen- 
dations to im.prove validation are provided. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 12 

A. BACKGROUND INFORMATION 12 

B. OBJECTIVES OF THIS WORK 16 

II. EXPERIMENTAL APPARATUS 1'^ 

A. INTRODUCTION 1'^ 

B. STEAM SYSTEM 1'7 

C. TEST CONDENSER 18 

D. CONDENSATE SYSTEM 18 

E. COOLING WATER SYSTEM 20 

F. SECONDARY SYSTEMS 20 

1. Vacuum System 20 

2. Desuperheater 21 

G. INSTRUMENTATION 21 

1. Flow Rates 21 

2. Pressure 22 

3. Temperature 22 

4. Data Collection and Display 23 

III. EXPERIMENTAL PROCEDURES 24 

A. OPERATING PROCEDURES 24 

1. Preparation of Condenser Tubes 24 

2. System Operation and Steady State 

Conditions 24 

B. DATA REDUCTION PROCEDURES 25 

1. Overall Heat Transfer Coefficient (U ) 25 

n 

2. Corrected Overall Heat Transfer 

Coefficient (U ) 25 

c 



5 



26 



3. Inside Heat Transfer Coefficient (h.) 

1 

4. Outside Heat Transfer Coefficient (h ) 26 

o 

C. DATA REDUCTION COMPUTER PROGRAM 27 

IV. RESULTS AND DISCUSSION 2 8 

V. CONCLUSIONS AND RECOMMENDATIONS 31 

VI. TABLES 32 

VII. FIGURES 39 

APPENDIX A: SAMPLE CALCULATIONS 57 

APPENDIX B: ERROR ANALYSIS 62 

BIBLIOGRAPHY 66 

INITIAL DISTRIBUTION LIST 67 



6 



LIST OF TABLES 



Table I. Location and Channels for Thermocouples 32 

Table II. Raw Data for Tube Number 1 33 

Table III. Tube Number 1 Results 33 

Table IV. Raw Data for Tube Number 2 34 

Table V. Tube Number 2 Results 34 

Table VI. Raw Data for Tube Number 3 35 

Table VII. Tube Number 3 Results 35 

Table VIII. Raw Data for Tube Number 4 36 

Table IX. Tube Niimber 4 Results 36 

Table X. Raw Data for Tube Number 5 37 

Table XI. Tube Nxxmber 5 Results 37 

Table XII. Data Points for Figure 14 38 



7 



LIST OF FIGURES 



Figure 1 . Droplet Path through a Tube Bundle with 

Side Drainage 39 

Figure 2. Photograph of Test Facility 40 

Figure 3. Schematic Diagram of Steam System 41 

Figure 4. Photograph of Test Condenser with Partial 

Insulation 42 

Figure 5. Test Condenser Schematic, Front View 43 

Figure 6. Schematic Side View of Test Tube 

Arrangement 44 

Figure 7. Schematic Diagram of Condensate 

and Feedwater System 45 

Figure 8. Schematic Diagram of Cooling Water System — 46 

Figiire 9. Schematic Diagram of Vacuum System 47 

Figure 10. Wilson Plot for Tube Number 1 4 8 

Figure 11. Wilson Plot for Tube Number 2 49 

Figure 12. Wilson Plot for Tube Number 3 50 

Figure 13. Wilson Plot for Tube Number 4 51 

Figure 14. Wilson Plot for Tube Number 5 52 

Figure 15. Comparison of Corrected Overall Heat 

Transfer Coefficient for Tubes 1, 2, 

3 with Fenner [5] 53 

Figure 16. Average Outside Heat Transfer Coefficient 

Ratio Versus No. of Tubes 54 

Figure 17. Sequence of Photographs Showing 

Side Drainage of Droplet 55 

Figure 18. Sequence of Photographs Showing 

Droplet Migration 56 



8 










JC-- 








NOMENCLATURE 



n 



C 

C. 



wv 



G 

GPM 

fg 

h. 

1 



n 



n 



Nu 

k 

w 

L 

LPM 

M 

ih 



m 

c 

n 

M 



cw 



2 

Outside Nominal Surface Area of Tube (m ) 

, 2 
Cross Sectional Area of Test Section (m ) 

Constant Determined from Experimental Data 

Specific Heat (kJ/kg-K) 

Cooling Water Velocity (m/sec) 

Inside Diameter (m) 

Outside Diameter (m) 

Drainage Factor 

2 

Mass Flow Rate of Cooling Water Per Unit Area (kg/m -sec) 
Gallons Per Minute of Cooling Water Flow 
Latent Heat of Vaporization (J/kg) 

2 

Inside Heat Transfer Coefficient (W/m -K) 

2 

Outside Heat Transfer Coefficient (W/m -K) 

th 2 

Mean Heat Transfer Coefficient on n^ Tube (W/m"^-K) 

Average Heat Transfer Coefficient over a Bundle 
of n Tubes 

2 

Nusselt's Heat Transfer Coefficient (V7/m -K) 

Thermal Conductivity of VJater (W/ra-K) 

Length (m) 

Liters Per Minute 
Wilson Plot Slope 

Mass Flow Rate of Condensate (kg/sec) 

Mass Flow Rate of Cooling Water (kg/sec) 

Number of Tubes 

Dynamic Viscosity (kg/m-sec) 



9 



Pr - Prandtl Number 

Q - Heat Transfer Rate (J/sec) 

Re - Reynolds Number 

3 

p - Density (kg/m ) 

R^ - Wall Resistance (m^-K/W) 

- Bulk Temperature (C) 

Tj^j, - Bulk Temperature (K) 

T . - Cooling Water Inlet Temperature (C) 

T^^ - Cooling Water Outlet Temperature (C) 

Tj^^ - Condensate Temperature in Primary Hotwell (C) 

"^shw ~ Condensate Temperature in Secondary Hotwell (C) 

T - Temperature of Steam Vapor (C) 

2 

- Corrected Overall Heat Transfer Coefficient (W/m -K) 

2 

- Overall Heat Transfer Coefficient (W/m -K) 

X - Abcissa of Wilson Plot 

Y - Ordinate of Wilson Plot (m^-K/W) 



10 



ACKNOWLE DGMENT S 



The author wishes to express his sincerest appreciation 
to his thesis advisor. Dr. Paul J. Marto, for his unfailing 
assistance and direction. A special note of thanks is 
deserved by Mr. Ken Mothersell for his technical and skill- 
ful support during this project. I would also like to thank 
Mr. Thomas Christian and LCDR Rick Foltz for their support 
and assistance. 

The entire faculty of the Mechanical Engineering Depart- 
ment is recognized for their loyal and untiring efforts in 
furthering my education to accomplish this goal. 

Finally, and certainly at the top of the list, I thank 
God for the talents and abilities he has entrusted to me and 
for the family he has surrounded me with, without whose loving 
support, sacrifices, and understanding this project could 
never have been completed. 



11 



INTRODUCTION 



I . 

A. BACKGROUND INFORMATION 

It has been well dociomented by Search [1] that naval 
marine condensers have been extremely conservative in design. 
Search concluded that heat transfer enhancement methods could 
decrease condenser space to weight ratio thereby establishing 
new design criteria for marine condensers. 

It has been the objective of past research by Beck [2] , 
Pence [3], Reilly [4], Fenner [5], Manvel [6] and Ciftci [7] 
to investigate the enhancement of heat transfer on single 
horizontal tubes. Their research concentrated on using differ- 
ent tube materials and various txibe geometries. Each of these 
investigations recommended that the effect of both inundation 
or condensate rain, and vapor velocity be examined for these 
enhanced tubes . 

Any attempt to understand the effect of inundation and 
vapor velocity must begin with the classical equation derived 
by Nusselt [8] for predicting the mean condensate heat trans- 
fer coefficient for a column of n horizontal tubes. It is 
normally written as a ratio of the average heat transfer 

coefficient for a column of n horizontal tubes (h ) to that 

n 

of the uppermost tube . The following relationship was 

obtained : 




n 



-1/4 



( 1 ) 
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where 



n = no. of tubes in a column. 

Experimental results to date indicate that Equation (1) 
underestimates the heat transfer. The above relationship 
was derived assuming; (1) condensation of a saturated vapor 
at negligible velocity and (2) laminar flow of the condensate 
film in a continuous sheet from one tube to the next at a 
constant temperature difference (as measured between the 
vapor and the wall) for all tubes in the bank. 

Neither of the above assumptions correspond to actual 
steam condenser conditions. Considering the first assump- 
tion, it is well known that with variable steam turbine speeds, 
the velocity of steam is often not negligible. Considering 
the second assumption, investigation, including filming, by 
Skklover and Buevich [9] reveals that condensate dropping 
onto the tube situated below in an "in-line" configuration 
does not drain as a continuous laminar sheet. In fact, it 
hardly spreads along the length of the tube at all, but rather 
rolls over it's perimeter forming droplets which locally 
thicken the film. Dr. David Eissenberg of the Oak Ridge 
National Laboratory in an unpublished v/ork preceding his 
Ph.D. dissertation [lOJ noted that the droplets do not only 
strike the tops of the tiobes below, but strike anywhere on 
the upper half. 

Eissenberg [10] performed experimental tests on a vertical 
column of five tubes in a triangular arrayed horizontal bank 
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with horizontal steam flow. His results were well above 
those predicted by Nusselt. He proposed a "side drainage" 
model with a droplet path as shown in Figure 1. His model 
assxjmes that some drainage occurs from the tube bottom to 
tube side in triangular arrayed tube banks thereby resulting 
in only part of the lower tube being inundated and the heat 
transfer coefficient for the tube being increased. If all 
the drainage was by the "side drainage" model, he showed 
that ; 



h 



n 



h 



= 0 . 6 + 0 . 42n 



-1/4 



Nu 



( 2 ) 



His experimental results fell between his side drainage 
prediction and that of Nusselt. 

A literature review reveals data scatter between that of 
Eissenberg and Nusselt. The data scatter results from the 
many experimental designs decided on by the individual authors. 
In summary, it can be generalized to say that in most cases 
the data has been fit to various modified forms of Equation 
(1) in which the exponent is empirically determined; 



h 

h 



n 

Nu 




(3) 



where values of s have been reported from 0.07 to 0.223. 

The effects of inundation and steam velocity are normally 
described in the literature as separate from each other. In 
fact, they occur simultaneously and their combined effect on 
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condensation in a tube bundle is complex. Fuji! [11] used 
the experiments of Nobbs and Mayhew to correlate the effect 
of inundation and vapor shear. The data for "in-line" tube 
banks resulted in Equation (4) for the first or uninundated 
tube : 



® 0 312 

Nu = 10.74 Re^^*-^ 

m L 



(4) 



where 



Nu = Nusselt number, h^d /k^ for pure 
m o o L 

steam without inundation; 

Re^ = two-phase Reynolds number, U d /v, ; 

L 00 O L 

k,. = thermal conductivity of liquid, 

1j 

d = outside diameter of tube; 
o 

V = kinematic viscosity of liquid, and 

JU 

is vapor velocity. 

O 

NUj^/Nu^ is computed for the data and inundated tubes resulting 
in Equation (5) : 



where 



O 



Nu /Nu 



m 



, , 0.071 (w/w°)°*^^ 

2 X 10° 



(5) 



Nu = mean Nusselt number for a tube 
m 

w = rate of inundation falling onto a tube 

w® = rate of condensation of a tube 

corresponding to Nu"* . 
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Fujii [11] determined that for the staggered tube bank, the 
inundation effect was smaller than that for the in-line bank. 

This would follow from examination of Eissenberg's side 
drainage model. 

Experimental data resulting from varying inundation rates 
and steam velocities is sufficiently scattered to suggest 
that no existing correlation adequately represents all the 
available data. This is because there are so many variables 
that affect condenser performance, some of which have yet to 
be explored (e.g., direction of steam flow, pitch to diameter 
ratio of the tubes, non-condensable gas effects, and the 
effects due to pitch and roll in a marine environment) . Research 
is presently being conducted in the United Kingdom aimed at 
producing a model of condensation suitable for sophisticated 
condenser performance calculations. However, the results have 
yet to be published. 

B. OBJECTIVES OF THIS WORK 

In order to evaluate the effect of inundation on a bank 
of enhanced tubes, it was necessary to establish data on a 
bank of smooth tubes as a standard of comparison. The objec- 
tives of this work were therefore: 

1. To modify the existing experimental apparatus to 
accommodate a bank of five tubes. 

2. To experimentally validate the system by comparing 
the experimental results to theoretical predictions, thereby 
establishing a data base on which to further study the various 
parameters affecting heat transfer within an enhanced tube bundle. 



16 



II. EXPERIMENTAL APPARATUS 



A. INTRODUCTION 

The existing test facility was designed by Beck [2] and 
built and tested by Pence [3] . In order to accomplish the 
goals of this work, major modifications were made to the 
original apparatus. In the following description of the 
experimental apparatus, attention will be focused on the 
changes and additions to the existing layout. 

B. STEAM SYSTEM 

The steam system is shown in Figure 3. The supply of 
steam is locally generated and supplied to the building which 
houses the experimental apparatus. The steam is provided via 
a 19.05mm diameter line and a steam inlet valve (MS-2). 

A pressure gage is located just prior to the steam separator 
which monitors the supply pressure as it is adjusted by 
(MS-2) . Following the steam separator, a line strainer pro- 
vides additional protection from contamination. After the 
strainer, the steam proceeds through a 31.75mm diameter line 
which provides for two possible steam paths. The primary 
path for system operation is via the throttling valve (MS-3) , 
through a desuperheater and into the test condenser. Inside 
the condenser, the steam is condensed on the test tubes. The 
steam which is not condensed proceeds via the vapor outlet 
on the test condenser to the secondary condenser. The secondary 
steam flow path is used to accomplish system stabilization 
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during start up and to control the mass flow rate of steam 
to the condenser during operation. Steam proceeds via 
(MS-4) directly to the secondary condenser. All steam 
lines except the primary path downstream of (MS-3) were 
insulated with 25.4mm thick fiberglass insulation. 

C. TEST CONDENSER 

The test condenser is shown from various views in Figures 
A, 5 and 6. Steam enters via the top and proceeds over the 
baffle separators and through a flow straightener , which is 
covered with three layers of 150 mesh screen, to the tube 
bundle. The condensate collects at the bottom of the condenser 
and flows out one of the two 12.7mm diameter openings at 
either end of the condenser to the hotwell where it can be 
collected and measxired. 

A considerable modification was made to the viewing sec- 
tion to accommodate observation of the condensation process. 
Three separate viewing windows each 203mm by 140mm by 12.7mm 
and made of pyrex plate glass were installed to provide 
maximum viewing of the active tubes while maintaining struc- 
tural integrity. 

The tube sheet arrangement is shown in Figure 6 . The 
tubes are arranged in a typical naval condenser spacing to 
diameter ratio of 1.5. There were five 15.9mm OD, 1.14mm 
thick, 304 stainless steel tubes that had cooling water 
passing through them. Although typical naval condenser tubes 
are made of 90-10 copper-nickel, the choice of 304 stainless 
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steel was based on "on-hand" stock and the fact that the 
principles of inundation do not rely on tube material although 
perhaps the heat flux may change due to a different tube 
wall resistance. The remaining half tubes were made of 
15.9mm OD aluminum bar stock and were fastened by screws 
to the outside wall of the steam flow path. This arrangement 
was selected to best simulate the steam flow conditions in a 
section of tubes within an actual condenser. The five test 
tTibes are singularly removable. The top tube can be replaced 
by a 304 stainless steel porous tube which could simulate 
various condensate inundation rates. 

The test condenser was insulated with a 25.4mm thick sheet 
of Armorflex insulation. 

D. CONDENSATE SYSTEM 

The . condensate system had no revisions and is identical 
to that shown in Figure 7 which is taken from Ref. [7]. The 
test condenser hotwell collects the condensate from the test 
tubes, while the secondary condenser hotwell collects the 
condensate from the secondary condenser. Valve (C-1) allows 
the isolation of the test condenser hotwell for condensate 
measurement. Since house steam was used as the steam supply 
system, the condensate collected in the hotwells is pxamped 
back to the house system by the condensate pump via valve 
(C-3) . The condensate lines were insulated with 19.1mm thick 
Johns-Manville Aerotube insulation. 
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E. COOLING WATER SYSTEM 



The cooling water system is shown in Figure 8. The 
water used was normal house water which had been passed 
through a water softener on the way to the supply tank. 

The water is pumped from the supply tank by a 5HP electric 
driven pump. It is routed to the flowmeter header via 51mm 
OD plastic pipe. The flow of cooling water for each test 
tube is then individually controlled by it's own rotameter. 
Each rotameter allowed a maximum flowrate of 70.4 LPM. The 
cooling water after passing through the test section was 
piped back to the supply tank. A separate system pumped the 
water through a filter and cooling tower returning the cooled 
water to the supply tank in an effort to maintain a constant 
cooling water inlet temperature. 

After leaving the rotameters, the system piping was 
reduced to 15.9mm, ensuring a distance of 1 meter ahead of 
the test section to ensure a fully developed velocity pro- 
file while passing through the test section. 

F. SECONDARY SYSTEMS 

1 . Vacuum System 

The vacuum in the test condenser and secondary con- 
denser was maintained by a mechanical vacuum pump and a 
vacuum regulator which induces air into the system. The 
vacuum piamp takes a suction from the secondary condenser hot- 
well which is connected to the test condenser hotwell via 
discharge piping. A cold trap at the inlet of the vacuum 



20 



pump forces incoming vapor to pass over a system of refrig- 
erated copper coils . This removes steam and entrained 
water from the vacuvim line preventing moisture contamination 
of the vacuum pump oil. The vacuum pump outlet is vented 
through a roof exhaust fan to avoid a health hazard from 
breathing any oil vapor exhausted by the pump. A schematic 
diagram of this system can be found in Figure 9 which is 
taken from Ref. [7] . 

2 . Desuperheater System 

The desuperheater removes sensible heat from the 
superheated system by injecting water at about 60 C via the 
existing feedwater system through valve (DS-1) and a rota- 
meter. The desuperheater is a 267mm diameter stainless steel 
can, 457mm high, having four nozzles inserted equidistant 
around the circumference of the inner top of the can. The 
nozzles are a fan type and are positioned such that the spray 
is downward at a 45° angle to allow for better mixing. A 
collection tank is located on the bottom of the desuperheater 
to allow for drainage of condensate. This system can be iso- 
lated by valve (DS-2) . 

G. INSTRUMENTATION 

1 . Flow Rate 

a. Cooling water flow rate was measured individually 
for the five separate tubes. Each flow rate was determined 
by a rotameter with a capacity of 70.4 LPM (18.6 GPM) . The 
calibration procedure used was identical to that listed in 
Appendix A of Ref. [3] . 



21 



b. Steam velocity was determined by calculation; 



Vel . 




(7) 



where 



m 



cond 



Mass flow rate of condensate 



A 



c 



(kg/sec) , = Q/h^^ 

2 

Cross sectional area (m ) , and 



V 



Specific volume of vapor (m^/kg) . 



2 . Pressure 

Two different pressure sensing devices were used 
during experimentation. They were a Bourdon tube pressure 
gauge which measured steam pressure and an absolute pressure 
transducer coupled with a 760mm mercury manometer which was 
used to measure test condenser pressure. 

3 . Temperature 

Stainless steel sheathed, copper-constantan thermo- 
couples were used as the primary temperature monitoring devices. 
Figure 5 shows the location of the five vapor thermocouples. 

The remaining 30 thermocouples of this type were located as 
shown in Figure 8, six on each tube, two measuring cooling 
water inlet temperature and four measuring cooling water 
outlet temperature. Calibration procedures of the thermo- 
couples were identical to those listed in Appendix A of 
Ref . [3] . 
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4 . Data Collection and Display 

An Autodata collection system was utilized to record 
and display the temperatures in degrees Celsius obtained from 
the primary stainless steel thermcouples and to record and 
display the pressure in cm of Hg inside the test condenser. 
Table I lists the channel numbers and location of these 
devices . 
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III. EXPERIMENTAL PROCEDURES 



A. OPERATING PROCEDURES 

1 . Preparation of Condenser Tubes 

Prior to installation, each tube was properly pre- 
pared to ensure filmwise condensation. The cleaning proce- 
dure is listed in Appendix A of Ref. [7] . 

2 . System Operation and Steady State Conditions 

The basic operating instructions developed by Pence 
[3] and modified by Reilly [4] were used. The only change to 
the procedures as listed in Appendix B of Ref. [7] was that 
instead of one cooling water flowmeter to adjust, the experi- 
menter had five to set as desired. 

In general it takes from two to three hours from 
initial light off until steady-state conditions are estab- 
lished. The parameters used in determining steady- state condi- 
tions were cooling water inlet temperature and steam inlet 
temperature. If the cooling water inlet temperature did not 
vary more than ± 0.6 C/HR or the steam temperature did not 
vary more than ± 0.3 C/MIN, state steady was considered 
achieved. 

The time for the system to stabilize between changes 
in cooling water flow rate during the Wilson Plot technique 
was approximately ten minutes. This time increment is suspect 
as other investigators waited about one hour for stabiliza- 
tion between changes. It must be pointed out however that 
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the amount of time required to collect data over five tubes 
in a system that cannot be shut down and repeated the next 
day prohibits the greater time increment between data points 
for the Wilson Plot. 

The general set up for the data taken in this research 
was a steam velocity of approximately 5m/sec, steam tempera- 
ture of 64 C and a test condenser pressure of approximately 
21 KPa. 



B. DATA REDUCTION PROCEDURES 

The raw data collected for each tube can be found in the 
Tables Section beginning on page 32* The following standard 
heat transfer equations were used to reduce the raw data into 
a form that can be used for evaluation. 

1 . Overall Heat Transfer Coefficient (U ^) 



U 



n 



1 " < 
n 



T 

sv 



T 



sv 




CO 



( 8 ) 



2. Corrected Overall Heat Transfer Coefficient (U ) 
c- 



U 

c 




(9) 



n 

where R^ is the wall resistance corresponding to different 
tube materials. This equation allows for the comparison of 
tubes of different materials for the same steam and cooling 
water condition within the test condenser. 
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3 . Inside Heat Transfer Coefficient (h j^) 

Nu = ^ = 0.036 Re°’^ (L/D)"°‘®^^ 

( 10 ) 

Equation (10) was selected because both the Dittus-Boelter 
and Sieder-Tate relationships which are commonly used assume 
a fully developed velocity, as well as, thermal profile. In 
this research, it was suspected that, although the velocity 
profile was believed to be fully developed, the thermal pro- 
file was not fully developed. When an L/D ratio of 57.6 is 
used in Equation (10) a constant of 0.029 results. This was 
validated by computing the average of the first three tube 
constants obtained as a result of the Wilson Plot technique. 
Wilson Plots for each tube are found as Figures 10 through 14. 

4 . Outside Heat Transfer Coefficient (h^ ) 



U 

n 



D 

p _ 

w D . h . 



( 11 ) 



The outside heat transfer coefficient is the parameter 
that is used to compare results of each tube in the bundle. 

Two very important assumptions were made in using this equation. 

a. The resistance due to fouling was equal to zero. 

This assumption is supported by the fact that the tubes were 
new, chemically cleaned and smooth. Also, treated soft water 
was used as the cooling medium. 
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b. The resistance due to non-condensable gases 
was equal to zero. This assumption was supported by the fact 
that the system was tested for air-tightness and found to be 
secure. In addition, it was believed that the velocity of 
steam passing through the test section was sufficiently large 
to keep the system purged of any non-condensables that might 
collect in the test section. 

C. DATA REDUCTION COMPUTER PROGRAM 

Reilly [4] developed the existing program in Fortran 
language. During this work, his program was translated into 
Basic language for use with the HP 9845 computer. This was 
done because of the author's desire to improve the speed and 
reliability of data reduction. Ultimately, this program with 
minor modifications, can be used in an integrated system 
between the Autodata Nine data collector and the HP 9845 com- 
puter. This will allow automatic data input with realtime 
data output for the experimenter. 
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IV. RESULTS AND DISCUSSION 



The validation of the experimental system was done first 
by comparing the first tiobe of this configuration with the 
smooth, 15.9mm OD tube of Fenner [5]. Secondly, the average 
heat transfer coefficient over n tubes divided by the heat 
transfer coefficient of the first tube for this experiment 
was plotted for comparison to the theoretical curves based 
on Equations (1) and (2) . 

The comparison of the number one tube of this research 
with that of a similar tube tested by Fenner [5] is shown in 
Figure 15. The discrepancies are attributed to the fact 
that the present system was operated at a steam velocity 
approximately three times the value used by Fenner, thereby 
causing the outside heat transfer coefficient to be higher 
during these tests. The general agreement v/as sufficient to 
validate the data input and reduction program. The plots 
for tube #2 and #3 show the expected decrease in overall 
heat transfer as one proceeds down through the bundle of 
tubes. This is caused by the decrease in vapor velocity as 
steam is condensed as it flows downward and the fact that the 
condensate film is thicker on the lower tubes due to inunda- 
tion from above. 

Figures 10 through 14 are the Wilson Plots that assist 
in determining the constant in Equation (10) . The first three 
tubes yield good linear plots with slopes which provide 
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constants of .026, .030 and .033 respectively. The slopes 
of the Wilson Plots for tubes #4 and #5 provided constants 
of .048 and .020 which were considered invalid. The data 
rrduction program gives the user the option of using the 
constant solved for via the Wilson Plot technique or inputing 
one of his own choosing. In this work, since the constant 
0.029 was the average over the first three tubes, it was 
used as input for all the tubes to determine the inside heat 
transfer coefficient which in turn was then used in the 
determination of the outside heat transfer coefficient. 

The ratio of h^/h^ as listed in Table XII was determined 
by taking the average outside heat transfer coefficient h^ 
for n tubes and averaging them, then dividing by the outside 
heat transfer coefficient of the first tube (h^^) . The results 
of Table XII are then plotted on Figure 16 along with the 
theoretical equations of Nusselt and Eissenberg. Based on 
the observations of condensate flow, shown in Figure 17 it 
was expected that the data over the first three tubes, where 
"side drainage" was the dominant mode of flow, would approach 
the Eissenberg curve. The data for tubes #4 and #5 was 
expected to fall closer to the Nusselt curve due to the observed 
presence of gravity dom.inated flow, but certainly not below 
it. 

The cause for lower outside heat transfer coefficients 
than expected may be due to several phenomena. During obser- 
vation of the condensate flow patterns on the first three 
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tubes, there was evidence of droplet migration, from the 
ends of the tubes towards the center. This migration is 
shown pictorially in Figure 18 and is presumably due to vapor 
flow axially along the tube. This vapor flow causes a non- 
uniform heat transfer rate across the length of the tube 
because of the increase in film thickness that develops toward 
the center of the tube instead of a thinner condensate film 
which would result if the vapor flow was uniformly downward. 
This thickening of the film in the central part of the tube 
would result in a lower average heat transfer coefficient 
than expected. 

The severe drop in performance of tubes number four and 
five indicate that either there is a very low mass flow of 
steam remaining to be condensed or there is a concentration 
of non-condensables that is affecting heat transfer performance. 
The latter is most suspect as the data for tubes four and 
five occurred between eight and twelve hours into the run 
allowing a significant buildup of non-condensables if steam 
velocity was not sufficient to purge it out. It is conceiva- 
ble that the non-condensables entered into the system from 
the steam side possibly from the house steam, the desuperheater 
spray system or via the vacuum bleed valve which was slightly 
open during operation to maintain desired test condenser 
pressure. The presence of secondary flow would also lead to 
a stagnation area which would trap non-condensables in pockets 
near the bottom and sides of the test section. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The experimental data obtained leads to the following 
conclusions : 

1. The system as modified is a valid test facility for 
further experimentation investigating the effect of innunda- 
tion and vapor velocity. 

2. There is evidence of secondary flow of steam within 
the test condenser which is suspected to have influenced the 
results . 

3. The measured ratio of is lower than expected. 

The following recommendations are provided: 

1. Improve the steam flow path to ensure a uniform 
downward profile through the bundle. This can be accom- 
plished by either a re-design of the steam inlet section or 
by a series of baffles spaced between the tubes to reduce 
the secondary flow effects. 

2. Measure the non-condensable gas concentration using 
the method described in Ref. [12]. 

3. Operate the system at higher steam velocities and 
mass flow rates to ensure adequate purging and heat transfer 
at the lower tubes. In addition, longer stabilization times 
between cooling water velocity changes should be tried in 
order to obtain more consistent Wilson Plots. 
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VI . TABLES 



Table I 

Location of Stainless Stell Sheathed 
Copper Constantan Thermocouples 
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Table II 



Raw Data for Ttibe Number 1 



% FLOW 


T . (°C) 

Cl 


T ^ (°C) 


O 

O 


GPM 




CO 


sv 


10 


22.5 


30.9 
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12.5 


22.9 
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Table III 
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Table IV 



Raw Data for Tiobe Niomber 2 



% FLOW 


T^. (°C) 

Cl 


T^^ (°C) 

CO 


T (°C) 

sv 


GPM 


10 


24.1 


31.8 


64.1 
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27.6 
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27.6 


64.2 
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Table V 





Results 


for Tube 


Nxomber 
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Table VI 



Raw Data for Tube Nxamber 3 



% FLOW 


T . (°C) 

cx 


T (°C) 

CO 


T (°C) 

sv 


GPM 


10 


24.9 
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Table VII 
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for T\ibe Number 
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Table VIII 



Raw Data for Tube Number 4 



% FLOW 
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T^^ (°C) 
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Table IX 
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for Tube Number 


4 




h. 


h 


C 


Re 


C 


1 


o 






wv 


6120. 


3149 . 


.029 


19662. 


1.17 


6904. 


2994. 


.029 


22820. 


1.36 


7660. 


3059. 


.029 


25945. 


1.55 


8427. 


3197. 


.029 


29249. 


1.75 


9125. 


3249. 


.029 


32228. 


1.94 


10510. 


3148. 


.029 


38351. 


2.33 


13168. 


2917. 


.029 


50706. 


3.11 


13168. 


2917. 


.029 


50706. 


3.11 


25131. 


2996. 


.029 


113553. 


7.00 


25131. 


2996. 


.029 


113553. 


7.00 



U 


U 


Q 


m 


n 


c 




c 


1789. 


1994. 


3054. 


.001300 


1804. 


2012. 


3118. 


.001328 


1882. 


2110. 


3259. 


.001387 


1985. 


2240. 


3437. 
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Table X 



% FLOW 



10 
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20.0 
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Raw Data for T\ibe N\Jir±)er 5 
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Table XI 
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for Tube 
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Table XII 



Data Results for Figure 14 



h /h 
n 



1 



n 



1 1 

.849 2 

.765 3 

.632 4 

.533 5 
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VII. FIGURES 




Fig. 1. Droplet Path Through a Tube Bundle with Side 
Drainage 
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Fig. 2. Photograph of Test Facility 





STEAM SYSTEM 




ro 

L 



41 



Schematic Diagram of Steam System. 
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Fig. 4. Photograph of Test Condenser with Partial Insulation 
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Test Condenser Schematic, Front View. 
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Schematic Diagram of Cooling Water System 



VACUUM SYSTEM 



I 



iA 

D 
< 
X 
OX 
K UJ 




U- H U 



u o: 
5 UJ 

cr O < 
IL I $ 



-P 

CO 

CO 



p 

o 

fO 

> 



o 

g 

03 

P 

CTk 

fd 

•H 

Q 

U 

•H 

4J 

td 



x: 

u 

CO 



<Ti 



CJ^ 

•H 



47 



I 



WILSON PLOT FOR TUBE #1 
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Fig, 10. Wilson Plot for Tube Number 



I 



M 



I 

1 



WILSON PLOT FOR TUBE *2 



O 

CO 



in 

< 

G 



CO 

CO 

(S3 

< 

CK 

CL 

OD 

d 

< 

Li 

(K 



I i_ 

CD 



(Sv 0 ix) Nfl/I 




04 



49 



Fig. 11. Wilson Plot for Tube Nxamber 
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Fig. 12. Wilson Plot for Tube Number 
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WILSON PLOT FOR TUBE #4 
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Fig. 13. Wilson Plot for Tube Number 
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WILSON PLOT FOR TUBE #5 
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Fig. 14. Wilson Plot for Tube Nuinber 
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hn^''h 1 




NUMBER or TUBES (n) 



Fig. 16. Average Outside Heat Transfer Coefficient 
Ratio Versus No. of Tubes 
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(c) (d) 

Fig. 17. Sequence of Photographs Showing Side Drainage 
of Droplet 
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Fig. 18. Sequence of Photographs Showing Droplet Migration 
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APPENDIX A 



Sample Calculations 

The following is an example of how the data reduction 
program progresses to the results. Tube number 1 at 25% flow 
was selected for this analysis. This same tube and flow rate 
was used for the error analysis in Appendix B. 

Input Parameters 

Outside Diameter (D ) 0.015875m 

o 



Inside Diameter (D^) 


0.0141m 


Tube Length (L) 


0.9144m 


Outside Nominal Surface Area (A^) 


0 .0456m^ 


Wall Resistance (Rw) 


5.72 X 10~^m^ - K/W 


Cooling Water Temperature Inlet 


24.1 C 


Cooling Water Temperature Outlet 


29.4 C 


Average Cooling Water Temperature 

(T, , T, , ) 

be bk 


26.75 C, 299.9 K 


Steam Vapor Temperature (Tg^) 


64.0 C 


Gallons Per Minute of Cooling Water 
(GPM) 


4.80 GPM 


Section 1. Water Properties 





y (MHUX’J) = (4.134 xlO”^) exp {[ (0 . 008291758) { 299 . 9) 



+ (2644.2189)/(299. 9) ] - 


10.59252566} 


-4 

y = 8.418 X 10 kg/m-sec 
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k(Kw) = 0.5565919 + (0 . 002174417) (26 . 75) - ( 0 . 7012 7 x lo"^) 

•(26.75)^ - (2.0914 X 10"^°) (26.75) ^ 
k = .6097355 W/m - C 

p (RHO) = 100.444434 - ( 0 . 1267336 8) ( 26 . 76) 

- (0.0023913147) (26.75)^ 
p = 999.343 kg/m^ 



C (CP) = 4.2377955- ( 0 . 0018553514 )( 26 . 75) 

+ (1.3948314 X lo"^) (26.75) ^ 

Cp = 4.198 kJ/kg - C 

m(MFRCW) = LPM x RHO x 1.67 x lO"^, where LPM = GMP x 3.78533 
= 18.169584 X 999.343 x 1.67 x lO"^ 

m = 0.30323 kg/sec 

Prandtl Number (Pr) 

Pr = uC /k = (8.418x10"'^ x 4 . 198 x 10^) /. 6097355 

P 

Pr = 5.796 
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Section 2 , Data Reduction 



1. Cooling Water Velocity (C ) = 4m/p7rD.^ 

WV 1 

= (4) (. 30323)7(999.343) (it) (.0141)^ 

1.9 43m/sec 



2. Mass Flow Rate Per Unit Area (G) 
G = = (p)(C^^) 

G = (999.343) (1.943) 

G = 1941.72 kg/m^ - sec 



3. Reynolds Number (Re) 

Re = G/u = ( .0141) (1941.72) 78.418 X 10 

Re = 32,523.46 



Determination of Overall Heat Transfer Coefficient (U^) 



U 



n 



m C T - T . 

E £n ^) 

A 'm — T 

n sv CO 



(.30323) (4.198 X 10-^) ,64 - 24.1, 

.0456 ^ ^64 - 29.4^ 



= 3978.63 W7m - C 
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5. Determination of Corrected Overall Heat Transfer 

Coefficient (U ) 
c 



U 

c 



U 



c 



U 



- R 



n 



w 



1 

3978.63 



_1 

- 5.72 X 10 



5150.85 W/m^ - C 



6 . 



Determination of Wilson Plot Parameters 
(a) Ordinate 

Y = J- = 1 

U 3978.63 

n 

Y = 25.13 X lo"^ m^ C/W 



(b) Abscissa 

Re°‘® Pr^/^ (32,523.46)°*®(5.796)^'^^ 

X = 13.67 X 10“^ 



7. Determination of Constant 



C = , where M = slope returned by linear regression 



subroutine 



M = .9735 



C = 



0.015875 



( .9735) ( .6097355) 



= .027 
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8. Determination of Inside Heat Transfer Coefficient (h.) 

1 

Nu = = 0.036 Re°*^ Pr^/^ (L/Do) 



where 



L/D = 



.9144 

015875 



= 57.6 



^i 



^ .029 Re°*^ Pr^/^ 



h. 

1 



029 (32,523.46)°'® (5.796)^^® ^'^^ ' 0141^ 



h. 

1 



= 9171.78 W/m - C 



9. Determination of Outside Heat Transfer Coefficient 



n. 



ife - - d7%- 



3978.63 



- 5. 72 X 10 



-5 



( .015875) 



(.0141) (9171.78) 



= 14,008 W/m - C 
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APPENDIX B 



Error Analysis 

The basic equations used in this section are reproduced 
from Reilly [4]. The general form of the Kline and McClin- 
tock [13] "second order" equation is used to compute the 
probable error in the results. For some resultant, R, which 
is a function of primary variables X^, the 

probable error in R, 6R is given by; 



12 n 

(B-1) 

where 6X, , 6X_., . . . , 6X is the probable error in each of 

12 n 

the measured variables. 

1 . Uncertainty in Overall Heat Transfer Coefficient, Un 



6U 6A „ 6C „ . . - 

^ = [(-^)^ + (-^)^+ (— )^ + ( 
'n ^n % m 



U. 



6T (T . - T ) 
SV Cl CO 



T -T . 

(T -T .)(T -T ) Jin 

sv Cl SV CO T___-T 



SV CO 



6T 



+ ( 



Cl 



T -T . 

(T -T .) <ln 
SV Cl T 



(■ 



6T 



CO 



-T 

SV CO 



T -T . 

SV Cl 

•sv ‘cd"“ T -T 

SV CO 



_)2jl/2 



(T -T Un 



The following values are assigned to the variables; 

2 

6A = ±0.0001 m 

n 

SC = ±.0042 KJ/Kg - C 

P 
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+.01 kg/sec 



6m 



6T 

6T 



sv 



Cl 



6T 



CO 



For Tube #1 at 



= ±.5 C 

= ±.l C 

= ±.l C 



25% Flow: 



6U 



n 



U, 



n 



r,^001,2 , ,.0042,2 , , .01 ,2 

^ .0456^ M.198^ \30323^ 



+ ( 



(.5) (-5.3) 



. (■ 



.1 



(39.9) (34.6) Jin 1.153' ^ M39.9)iln 1.153 



r) 



+ ( 



) 2 ] 1/2 



34.6 Jin 1.153 



6U 



n 



U, 



= .044 



n 



Uj^ Tube #1 - 25% = 3978 ± 175 W/m^ - C 



2 . Uncertainty in Inside Heat Transfer Coefficient (hi) 



The probable error in the inside heat transfer coefficient 
is given by: 



^ ^i r,6k,2. /’\,2 . ,0.86Re,2. ,0.3336Pr,2. ,60,2^1/2 

"h~ = + (-dT^ ( Re' ^ ^ ^ ^ 

where 



6k = 


±0.001 


K/m - C 


6Di = 


±0.001 


m 
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6Pr = ±0.10 

6C = ±0.001 



6 Re 
Re 



[ (^) ^ + 






i 



[(. 01 )^ + ( 



10 



8.418 



2 



+ 



,.001 , 2 , 1/2 

\0141^ J 



= .073 



&h. 
1 






.6097' 



.0141' 



.073) 2 + ( -333 .1, 



5.796 



. 001 . 2 1/2 

\ 029 ^ ^ 



= .098 



hj_. Tube #1 25% = 9171 ± 901 W/m^ - C 



3. Uncertainty in the Outside Heat Transfer Coefficient ( h^) 

The probable error in the outside heat transfer coefficient 
is given by: 



h 



6U 



= { [- 



n 



-]^ + [- 



6Rw 



- Rw - 5^) 

D 6h. 

(5-^) (17:'-) 

^ - D 

n 11 



n 



- Rw - 



D : h . ' 
1 1 
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where 



<5U 



n 



U 



044 



n 



6Rw = 2.86 X lo"® - C/W 



<5h. 
1 

hi 



= .098 



1 T,.. Do 

— — RW — -zrT — T— r- 
Dl hi 



7. 14 X 10~^ - C/W 



ho 



= { [- 



044 



(3978.63) (7.14 x 10 



,]2 + ^ 2. 86 X 10~^ j 2 



7.14 X 10 



-5- 



6h 



^ ^ (1.2276 X 10~^) ( . 09 8) j 2j 1/2 
7.14 X lO"^ 



= .174 



h^ = 14,008 ± 2435 W/m - C 
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